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Abstract 

The inverse ac Josephson effect involves rf-induced (Shapiro) steps that cross over the zero-current 
axis; the phenomenon is of interest in voltage standard applications. The standard analysis of the 
step height in current, which yields the well-known Bessel- function dependence on an effective ac 
drive amplitude, is valid only when the drive frequency is large compared with the junction plasma 
frequency or the drive amplitude is large compared with the zero- voltage Josephson current. Using 
a first-order Krylov-Bogoliubov power-balance approach we derive an expression for the threshold 
value of the drive amplitude for zero-crossing steps that is not limited to the large frequency or 
large amplitude region. Comparison with numerical solutions of the RSJ differential equation shows 
excellent agreement for both fundamental and subharmonic steps. The power-balance value for the 
threshold converges to the Bessel- function value in the high-frequency limit. 
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The term 'inverse ac Josephson effect' was coined by Levinsen et al. jD to describe Shapiro 
steps 0] in the current-voltage characteristic of a small Josephson junction that cross over the zero- 
current axis. The phenomenon was recognized as being potentially important for voltage-standard 
applications inasmuch as the elimination of the dc bias current would eliminate one possible pathway 
for the entry of noise into the system. It was studied extensively — experimentally analytically 
and computationally — by a number of authors (see, e.g., Kautz || plus references therein). The 
'standard' analysis that emerged from these studies expressed the height in current of phase-locked 
Shapiro steps in terms of Bessel functions of an effective drive amplitude; this analysis was shown 
to be appropriate when the frequency of the ac drive is large compared with the junction plasma 
frequency or when the amplitude of the ac drive is large compared with the zero- voltage Josephson 
current. Our analysis, instead, is based on the power-balance formalism developed by Krylov and 
Bogoliubov Q; it is not limited to the large frequency or the large amplitude regions. Comparison of 
the predictions of the power-balance approach with results obtained by direct numerical integration 
of the model differential equation typically show agreement to at least three significant digits. In 
the high-frequency limit they approach those predicted by the Bessel-function expressions. 

The starting point of our analysis is the usual shuntcd-junction model of the small Josephson 
junction, subjected to a dc-bias current and an ac-driving current; the corresponding differential 
equation, in normalized form, is 

4> + c«j) + sin <f) = r + 7 cos(uit) , (1) 

where T and 7 are, respectively, the dc- and ac-driving currents. 

The corresponding unperturbed equation, i.e., with a, T, and 7 set to zero, has two types of 
solutions, viz., oscillatory and rotary. Since we are interested in non-zero-voltage states of the 
junction we focus on the latter, which have the form 

<K*)=2am(^;k), (2) 

where am() is the Jacobian elliptic amplitude function of modulus k, and t Q is an arbitrary 
constant. The instantaneous junction voltage in this case is thus given by 

^) = ~dn(^;k). (3) 

We now return to the original model given by Eq. (1). An obvious physical requirement for 
the junction voltage to be a stationary, periodic function similar to that given by Eq. (3) is that 
the average power dissipated, i.e., the conductance times the mean-square voltage, be equal to the 
average power furnished by the drive currents, i.e., the mean- value of the drive current times the 
voltage. Assuming that Eq. (2) can be used as a first approximation to the solution of Eq. (1) in 
the presence of dissipation and drive, we can write, using results from 0], the time-average power 
dissipated as 

I2 AaE(k) 

p - = a< ^ >= Fi# (4) 

where K{k) and E{k) are, respectively, the complete elliptic integrals of first and second kinds. 
Assuming first the presence of only a dc-bias current, i.e., setting 7 = 0, we can write the time- 
average input power as 

7rr 

P^o = r<d»= m - y (5) 
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Equating Eq. (4) to Eq. (5) gives an expression for the McCumber branch of the current-voltage 
characteristic of the junction in the following parametric form 

AaE(k) 

Vs< +>=kwy ^ 

In the zero-bias configuration, instead, i.e., with T = and 7 7^ 0, we can write the input power as 

Pin,ac = 7 < 4> COs(wi) > , (8) 

which, by expressing the Jacobian dn() function in terms of its Fourier-series expansion we can 
write as 

_ 7T7COs(a;f) 4tt7CQs(W) ^ q m mirjt - t Q ) 

t\ n ,ac-< kR{k) + kR{k) ^ 1 + q 2m C08 i kK{k ) J>> W 

where q = exp(—irK(k')/K(k)), in which k' is the complementary modulus j^]. From Eq. (9) we 
see that Pi n ,ac 7^ only if 

"=wm (10) 

for some integer m, and, assuming Eq. (10) to be satisfied, that Pin.ac can vary smoothly from 
zero to (Pin,ac)max depending on the value of the phase-shift term t a (we have P in ,ac = (Pin,ac)max 
when t Q = 0). The threshold value, Jthr, of the drive amplitude is the minimum value for which it 
is possible to satisfy the equation P out = {Pi n ,ac)max', this yields the value 

2aE(k) 1 + q 2m 

lthr = ; — • (11) 

irk q m 



For 7 > jthn the constant t adjusts itself according to the relation 
, mirt ~f thr 



(12) 



Finally, in the general case in which both T 7^ and 7 ^ 0, the dissipation defined by Eq. (4) is 
balanced by both dc and ac inputs. This permits, e.g., calculating the minimum bias-current value 
to which a phase-locked step extends: assuming Eq. (10) to be satisfied, we find 

_ 4aE(k) 2 iq m 

mm ~ irk l + <? 2 ™' 1 ' 

which, for 7 > j t hr [see Eq. (11)], is a negative number. 

The 'standard' Bessel function expression for the height in current of a phase-locked Shapiro 
step is H § 

Ar = 2J„(^^=), (14) 
wyw + or 

where J n Q is the Bessel function of order n. Thus, from Eqs. (6) and (13), the threshold value of 
the drive amplitude in this approach is 

u ] thr - B ) = ^m. ( i5) 

ojVoj 2 + a 2 nk 
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At this point it must be mentioned that the integer n in Eq. (14) is not the same as the integer 
m in Eq. (10): in fact, the Bessel function index n refers to a superharmonic number whereas the 
power balance index to refers to a sM&harmonic number. To be perfectly clear on this point, the 
corresponding Shapiro steps in the current- voltage characteristic of the junction occur at normalized 
voltages 

77 

V =<cf»=— w. (16) 

TO 

Thus, the only place where the two expressions, Eq. (11) and Eq. (15), can be compared directly is 
at the fundamental frequency, where n = to = 1. 

Our results are summarized in Fig. la,b. In these figures, the solid curve is the power-balance 
threshold, calculated from Eq. (11), the dashed curve is the Bessel-function threshold, calculated 
from Eq. (15), and the diamonds are values obtained from the direct numerical integration of Eq. 
(1). Except for the highest point in Fig. lb, where the numerical result begins to diverge slightly 
from the power-balance prediction, these two typically agree to at least three significant digits. A 
similar agreement between power-balance prediction and numerical experiment is obtained at the 
second subharmonic (to = 2). For higher frequencies than those shown in Fig. la, the two predicted 
threshold values, Eq. (11) and Eq. (15), converge to a common asymptotic limit. 

In order to calculate numerically the points indicated by diamonds in Fig. 1, we often found it 
necessary to 'tune' fairly precisely the initial conditions used for Eq. (1). This fact suggests that, at 
least in some cases, the basin of attraction of the phase-locked state is rather small, perhaps indeed 
vanishing, a fact that would have important implications for the stability of the locked state, e.g., 
against thermal fluctuations. The technique of cell mapping |7j] might be a useful tool for exploring 
this question. 

Braiman et al. have extended the Bessel function approach to subharmonic frequencies. In 
fact, their Eq. (2) has just the form of the Fourier-series expansion of the Jacobian am() function 
of our Eq. (2) §. This approach gives expressions for the step height in current — corresponding to 
Eq. (14) for superharmonic steps — in terms of sums of products of Bessel functions. The approach, 
however, clearly becomes rather unwieldy if more than the first few terms in the Fourier expansion 
are employed. 

The power-balance approach can also be extended beyond the lowest level of approximation that 
we have employed here J^]; this would presumably permit a description of superharmonic steps, 
which presently is lacking from our analysis. However, superharmonic steps occur at progressively 
higher voltages, which implies progressively higher frequencies, where the simple Bessel-function 
expression, Eq. (15), is known to give a reasonable description of the situation. Thus, the reward 
to be obtained here might not be worth the effort. 

Finally, we must offer one caveat: as is apparent from Fig. la, the power-balance prediction 
becomes progressively better than the Bessel-function prediction as the ac drive frequency is reduced 
below unity (which is the plasma frequency, with our normalization). This, however, is a region 
where, in addition to simple periodic solutions of Eq. (1), exemplified by the ansatz of Eq. (2), 
there are known to exist also complicated quasi-periodic and chaotic trajectories |9j. Consequently, 
whereas Eq. (11) does give a good estimate of the threshold value of the drive amplitude if a simple 
step exists, it does not guarantee the existence of such a step. 
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Figure Captions 

Fig. 1. Comparison of power-balance prediction (solid curve) and Bessel- function prediction (dashed 
curve) with numerical integration result (diamonds), at fundamental frequency (n = m = 1). 
Parameters: (a) a — 0.01; (b) uj — 1.0. 
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